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Abstract. By developing the Tanaka theory for rank 2 distri- 
butions, we completely classify classical Monge equations having 
maximal finite-dimensional symmetry algebras with fixed (albeit 
arbitrary) pair of its orders. Investigation of the corresponding 
Tanaka algebras leads to a new Lie-Backlund theorem. We prove 
that all flat Monge equations are successive integrable extensions 
of the Hilbert-Cartan equation. Many new examples are provided. 



Introduction and main results 

In his paper of 1895 [L], S. Lie made the important observation that 
the simplest differential equations, beyond those of finite type, are 
equations with general solution depending upon 1 function of 1 vari- 
able, and he indicates that they are integrable by ODE methods. 

Such differential equations can generally be encoded as rank 2 dis- 
tributions or equivalently as codimension 2 Pfaffian systems on k- 
dimensional manifolds. For = 5, Goursat showed [G] that these 
distributions can be realized as Monge equations 

y' = Fix,y,z,z',z"), (1) 

while E. Cartan jClj realized them as reductions of involutive pairs of 
2nd order PDEs in 1 dependent and 2 independent variables. The inter- 
play between the Monge and PDF realizations of rank 2 distributions 
is an important aspect of their geometry. 

From the Goursat-Monge viewpoint, one can therefore interpret Car- 
tan's 1910 paper as a detailed analysis of the geometry of Monge equa- 
tions ([1]). In what is still considered to be a computational tour de 
force, Cartan solved the equivalence problem for ([1]) and, as a corol- 
lary, deduced that the celebratecOl Hilbert-Cartan equation 



y' = {z"Y (2) 



^We write this in retrospection [HllC2| : in |C1| only the PDE model was explicitly 
constructed. 
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is the Monge equation of type ([I]) with the largest finite-dimensional 
Lie algebra of symmetries. 

Despite Lie's observation and Goursat's contribution, or perhaps be- 
cause of the complexity of Cartan's analysis, little progress has been 
made in the symmetry analysis of the general Monge equations 

l/("^) = F(x,i/,i/',...,y("^-^),^, (3) 

(here and beyond y = y{x), z = z{x) and we assume m <n). 

These underdetermined ordinary differential equations are important 
in the variational calculus and control theory, in reductions of compat- 
ible PDEs and in the symmetry approach to Darboux integrability. 
Their geometry is the subject of the present article. 

The first goal of the paper is to describe the most symmetric equa- 
tions of type ([3]). There are different kinds of symmetries and Lie- 
Backlund type theorems clarify the relations between point, contact, 
generalized (higher) and internal symmetries whenever possible, see 

[un iKLviKKOick] . 

Below we restrict to purely internal symmetries. We will assume 
m > 0, otherwise ([3D has an infinite-dimensional algebra of symmetries. 

Theorem 1. Among equations £ of type ^ with fixed m < n and 
F^{„)^{n) 7^ 0, the most symmetric (realizing maxdimsym(£^)j are in- 
ternally equivalent to the equation 

The dimension of the symmetry algebra sym(£^) for this equation equals 
2n + 5 for m = 1, n > 2 and 2n + 4 for m > 1. 

Let us denote by £m,n the model equation (jlj). Notice that the equa- 
tion £^1^1 is internally equivalent to the pure jets J^(R, R) with (infinite- 
dimensional) contact symmetry algebra, while Si^2 is the Hilbert-Cartan 
equation ([2]) with symmetry algebra 02 = Lie((j2) of dimension 14. 
These are the exceptional cases omitted in Theorem [TJ 

Remark 1. As a by-product of the proof of TheoremUl we conclude 
that Monge equations £ with F^(n)^(n) ^ are strongly regular in the 
sense of Tanaka [Tj. This implies the existence of a Cartan connection 
for ^ as a im,n-valued 1-form on a natural bundle over the equation 
£, where tm,n the Tanaka algebra calculated in Section ISTB . 

As is well-known, equation ([3]) defines a submanifold £ in the space of 
jets J'"'"(Ri, ^ M'»+'^+3(a;^ y,yi, . . . ,y^, z, zi, . . . , z^), with dim£ = 
m + n + 2. The canonical rank 3 distribution 

C = Ann{(iyj — yj+i dx, dzj — zj^i dx\l<i<m,l<j<n} 
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(Cartan distribution) on the jet space J'"'"(]R^, M^) induces a rank 2 
distribution A on Internal symmetries of are, by definition, 
symmetries of the pair {S,A). 

As noticed above, rank 2 distributions were investigated in many 
classical papers. Recently Doubrov and Zelenko pZj described their 
internal geometry via a canonical frame on a bundle over the manifold 
S, providing a unique solution for maximally symmetric model among 
all rank 2 distributions from some generic stratum (called maximal 
class in loc.cit). This class does not appear to have a clear description 
for distributions, arising for differential equations 

Theorem [T] yields maximally symmetric distributions in a large class 
of distributions corresponding to Monge equations of a fixed dimension. 

Corollary 1. Among all rank 2 distributions A corresponding to equa- 
tions ^ with fixed m + n = k + 1 and F^(n)^{n) ^ 0, the most symmetric 
are locally equivalent to the canonical distribution of the equation Si^k- 

The non-degeneracy condition F^(n)z(n) 7^ is related to the finite- 
dimensionality of sym(£^) and will be explained in Section [21 It im- 
plies that the growth vector of the corresponding distribution starts 
(2, 1, 2, ... ) and that the distribution is not de-pro longable. 

Let us note that, according to |AKOj . internal symmetries for equa- 
tions of type ([3]), in the particular case m = n = 2, coincide with 
first order generalized symmetries. We can prove much more than this, 
namely, we prove a new Lie-Backlund type theorem. 

Theorem 2. For equations of type with F^(n)^(n) ^ 0, < m < n 
and {m,n) being different from the exceptional cases (1,1) and (1,2), 
the algebra of internal symmetries coincides with the algebra of external 
symmetries, induced by Lie transformations of the ambient mixed jet 
space J'^''^(M,M2). 

Lie transformations of J''"'"(M, M^) will be described in Section |H Let 
us remark that Theorem [2] is the first instance of Lie-Backlund type 
theorem whose proof entirely depends on the structure of the graded 
Lie algebra associated to the distribution, and cannot be derived by 
the usual de-prolongation arguments. 

To prove our results we exploit and elaborate upon the Tanaka the- 
ory [T]. Distributions corresponding to equation (jl]) are fiat in its 
sense. There exist however underdetermined ODE systems, for which 
the respective rank 2 distributions cannot be realized as non-degenerate 
Monge equations. Their fiat models can nevertheless be calculated via 
the same technique. 
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In Section [5] we relate the algebraic theory of central extensions of Lie 
algebras to the geometric theory of integrable extensions. We provide 
a very simple and constructive description of all flat models. 

Theorem 3. Any flat rank 2 distribution with growth vector (2, 1,2,...) 

(be it represented as a nondegenerate single Monge equation or underde- 
termined system of ODEs) is a result of successive integrable extensions 
of the Hilbert-Cartan equation 

This allows us to calculate many examples of flat rank 2 distributions 
and the corresponding ODE systems. 

At the end of the paper we show that our theory does not restrict 
to ODEs only. We demonstrate relations of our results to integration 
of determined and overdetermined PDEs. In particular, we construct 
a flat non-trivial integrable extension of Cartan's most symmetric in- 
volutive PDE model of 1910. 

The results of this paper indicate that the method of E. Cartan from 
his seminal paper |C1] can be extended to a much broader class of PDE 
systems. This generalization and further applications to the method of 
Darboux [AFV] and compatible overdetermined systems |KL] will be 
presented in a separate publication. 

Acknowledgment. We are grateful to Valentin Lychagin for dis- 
cussions during the initial stage of the project. Hospitality of Institut 
Mittag-Leffler (Stockholm) in 2007, MFO (Oberwolfach) in 2008 and 
Banach center (Warsaw) in 2009 is acknowledged. LA. was supported 
by NSF grant DMS-0713830. The computations in Sections [5] and [6] 
were done, in part, using the Differentia I Geometry package of Maple. 

1. Tanaka algebra of a distribution 

Given a distribution A C TM, its weak derived flag {Aj} is given via 
the module of its sections by L(Aj+i) = [L(A),L(Aj)] with Ai = A. 
The strong derived flag is given by the filtration Vj+i = dVi with 
F(V,+i) = [L(V,),L(V,)] and Vi = A. 

We will assume throughout this paper that the distribution is com- 
pletely non-holonomic, i.e. A^ = TM for some k, and we also assume 
that A is regular, so that the ranks of both Aj and Vj are constant. 

The quotient sheaf m = ©t<o5'i, 9% = A_j/A_i_i has a natural struc- 
ture of the graded nilpotent Lie algebra (GNLA) at any point x G M. 
The bracket on m is induced by the commutator of vector fields on 
M. A is called strongly regular if the GNLA m(x) does not depend on 
x G M (but we do not impose this property a priori). 
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The growth vector of A = g^i is the sequence of dimensions 
(dim^f.i, dim(7_2 . . .)• 0-"^ other sources, the growth vector is one of 
the sequences (dim Ai, dim A2, . . . ) or (dim Vi, dim V2, • • • ) but, as our 
main tool is the Tanaka theory, we choose this convention). 

The Tanaka prolongation g = m is the GNLA with negative graded 
part m, whose non-negative part is defined successively by the rule [Xl 
lY] : gk consists of all m G gk+i^g*, acting as derivations: u{[X, Y]) = 

i<0 

[m(X), Y] + [X, u{Y)], X,Y e m. Notice that since m is fundamental, 
i.e. g^i generates the whole GNLA m, the grading k homomorphism u 
is uniquely determined by the restriction u : g^i gk-i- 

The space g = (B gi is naturally a graded Lie algebra, which we will 
call the Tanaka algebra of A (it depends on the point x G M, unless 
A is strongly regular). By contrast, the GNLA m = ©i<ofi'i will be 
called the Carnot algebra of A (it is also called Levi- Tanaka or symbol 
algebra in the literature). 

We wish to indicate that g = rh can alternatively be defined via 
right extensions for modules, just as right extensionqj for Lie algebras 
are governed by H^{m, m) [SLt lF]. First notice that the cohomology W 
of GNLA is graded. We denote the grading by a subscript W = ®H]^. 

Then g^ = ifQ(m, m) and we define g° = m © go. On the next step 
we get gi = iJ^(m, 0°) and we define fl^ = 0° © gi. This is equivalent 
to adding variables Wi to kill the cohomology classes [ui] G if^(m, g°), 
that is, let dwi = Ui. The m-action naturally extends to making 
it into a representation and we continue: 0^ = g^k ® ■ ■ ■ ® gp with 
gp = Hp{m, g^'^). The prolongation of m is g = limg^. 

The importance of the Tanaka algebra is in majorizing the symmetry 
algebra sym(A) of the distribution. The following is the Corollary of 
Theorem 8.4 from [T]: 

Theorem 4. Let (M, A) be a manifold equipped with distribution, m 
its Carnot GNLA and g = m the Tanaka algebra. Then 

dimsym(A) < dimg. 

Moreover, provided g is finite- dimensional, this inequality is an equality 
if and only if the distribution A is locally flat. 

Some remarks concerning this theorem are of order. First, this theo- 
rem is formulated (and proved in [T]) for strongly regular systems, and 



Both fundamental extension constructions in the theory of Lie algebras play a 
role in this paper: right extensions are important for Tanaka prolongation, while 
central extensions appear in Section [5] in connection with integral extensions. 
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this will suffice for our purpose. But it remains true if we understand 
the right hand size as sup^ dim0(a;). 

Second, the last part of the theorem is not formulated in [T] as we 
have, but it follows from the construction of an absolute parallelism on 
the prolongation manifold of the structure. 

Finally, locally flat means that the structure functions of the absolute 
parallelism vanish. This is equivalent to the claim that the distribu- 
tion is locally diffeomorphic to the standard distribution on the affine 
realization of Lie group corresponding to the GNLA g, see 

2. Finite-dimensionality of the symmetry algebra 

In this section we give a sufficient condition for dimg, and hence 
dimsym(A), to be finite in the case dim A = 2. The statement has 
been surprisingly unknown until now. 

Notice that for a rank 2 non-holonomic distribution the growth vector 
always begins like (2, 1, ... ) and that the ffist number after a sequence 
of I's can be only 2. 

Theorem 5. Let A be a completely non-holonomic rank 2 distribu- 
tion with growth vector (2, 1, . . . , 1, 2, . . . ). Then the algebra sym(A) is 
finite- dimensional. 

Proof. Let g = (Bgi be the Tanaka algebra of A evaluated at any point 
X G M. In view of Theorem [Hit is enough to verify that dimg < oo. 

To this end we use Corollary 2 of Theorem 11.1 from [TJ, which 
states that if the subalgebra 

i)o = {v ego: [v, gr] = Vr < -1} 

is zero, then q is finite-dimensionalU 

Let —k be the grading of the second two-dimensional space g-k in 
m (fc > 3). Choose a basis ei,e[ of g-i that generates basis elements 
62 = [ei, e[] of g^2, es = [ei, 62] of g^s, Ck-i = [ei, ek-2] of gi-k and 
Cfc = [ei,efc_i], e'fc = [e[,ek-i] of g^k- 

Let f G P)o and write [v,ei] = aci + be[, [v,e[] = cei -\- de[. Then 
[v,efc_i] = implies [v,ek] = [aci -\- be[,ek-i] = ack + be'^ = and 
[v, ey = [cei + de[, Ck-i] = cck + de'i^ = 0. Whence a = b = c = d = 
and therefore v = 0. The claim follows. □ 

Consider the growth vector of A. If it starts (2, 1,2,...) then, by 
Theorem 0, the symmetry algebra is finite-dimensional. The other 



^The meaning of the condition f)o = is this: the complex characteristic variety 
of the Lie equation for sym(A) is empty. 
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possibility, namely (2, 1, 1, ... ), has a clear geometric meaning which 
we provide in the next theorem. 

Recall that the (geometric) prolongation of a rank 2 distribution A 
on M is (we consider only the local picture) the rank 2 distribution 
A = (cos t-U + sin t-V,dt) on M = M x , where U, V are generators 
of A and t the coordinate on S^. 

If the distribution A has preferred (vertical) section V G r(A), then 
the (affine) prolongation A = A*^^) is simply A = {U + tV,dt) on 
M = M X M^, where t is the coordinate on R^. This coincides with 
(Spencer) prolongation in the geometric theory of PDEs ( [KLVt [MZ] ) . 

The reverse operation is not generally possible but, in the case when 
A = A^^\ we shall call the operation A t— A de-prolongation. 

Theorem 6. Let A be a rank 2 non-holonomic distribution with the 
derived rank 3 distribution A2 = [A, A]. The growth vector of A is 
(2, 1,1,...) if and only if there exists a vector field ^ G r(A) which is 
a Cauchy characteristic for the distribution A2. In this, and only in 
this case, A can be locally de-prolonged. 

Recall that the Cauchy characteristic ^ is a symmetry of the distribu- 
tion A and its derived flag. So de-prolongation in this case is simply the 
quotient of the first derived distribution by the Cauchy characteristic 

(M,A) = (Af, A2)/e 

The above theorem is due to E. Cartan [C2] . We present a new proof 
here for the convenience of the reader. 

Proof. The claim about the existence of a Cauchy characteristic is ob- 
vious. We shall only demonstrate that A*^^) = A. 

For this, let us write M locally as the product Wj:^^ x in such a 
way that ^ = df. We take the other generator X of A from TR^~^ with 
coefficient dependent on {x,t). The derived flag equals 

A, = A={X,dt), A2 = {X,X,dt), A, = {X,X,X,[X,X],dt), 

where dot denotes the differentiation by t (or equivalently, the Lie 
derivative along dt). Since the growth vector is (2, 1,1,...) and dt is 
the Cauchy characteristic, we must have 

X = aX + (3X 

for some t-dependent functions a,(3 e C°°(M x M). Let W = u X + v X 
with M, f G C°°{M X R). Then solution to the system 



W = {ii + av)X + {u + i! + f3 v)X = 
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is equivalent (provided rank(X, X) = 2) to the ODE system 

dt \v) - \-l -P) \v 

By choosing a basis of fundamental solutions, we get two t-independent 
vector fields U,V such that {U,V) = {X,X). Since generators are 
defined up to scaling, we can assume X = U + bV . Regularity implies 
b ^ 0. By reparametrizing t we can assume that b = t. Thus we have: 

A={U + tV,dt), A2 = {U,V,dt), A = {U,V), 
so that the claim follows from the definition of the prolongation. □ 

De-prolongation does not just result in the removal of 1 from the 
growth vector at the first component and shift. It can significantly 
change the weak derived flag. 

Example. Consider the Hilbert-Cartan equation £^12. It has growth 
vector (2, 1, 2) and the Tanaka algebra is the exceptional Lie algebra Q2 
with grading (the number in the sequence below means the dimension 
while the subscript - the grading): 

(2(-3), l(-2),2(-i),4(o),2(i), 1(2), 2(3)). 

The prolongation of £^1,2 is a 6-dimensional manifold with rank 2 dis- 
tribution of growth (2, 1, 1, 1, 1). Here we see that de-prolongation is 
not just removal of 1 and shift. 

The Tanaka algebra of this 6-dimensional model is still 02, but now 
with the grading 

(l(-5), l(-4), l(-3), l(-2), 2(-i), 2(0), 2(1), 1(2), 1(3), 1(4), 1(5))- 

The next prolongation has growth vector (2,1,1,1,1,1), and at this 
point the Tanaka prolongation is the contact (infinite-dimensional) al- 
gebra. This means that the 2nd (Spencer) prolongation of Si^2 is not 
Tanaka- fiat. 

Since the algebras sym(A) and sym(A*^^'') are isomorphic. Theorems 
[6] and [5] imply the following important statement (it can be also deduced 
in a non-trivial way from Theorem 7.1 of |MTj ). 

Corollary 2. A regular germ of a rank 2 distribution A has infinite- 
dimensional symmetry algebra only in one of the following cases. 

(1) A completely non-holonomic: then it is equivalent to the Cartan 
distribution Ck on jet-space J'^(R, R), with k = dim M — 2. 

(2) Aj 7^ TM for any i: then there exists a codimension 1 foliation 

of M with TT D A such that all leaves (JF, A) are equivalent. 
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Recall that distributions A with strong growth vector (2, 1, 1, . . . , 1) 
are called Goursat distributions [G]. Their regular points can be char- 
acterized by the condition that the growth vector is (2, 1, 1, . . . , 1) in 
both the weak and the strong sense. Near such points the normal form 

Ck = Ann{dyi - y,+i dx\^ < i < k} C TJ''{R^, M^) 

is provided by the von Weber - Cartan theorem [WllC2j . see also |GKR] . 

This concludes our discussion of distributions with growth vector 
(2,1,1,...). From this point forward we assume that the distribution 
is not de-pro longable, i.e. its growth is (2, 1, 2, ... ). 

3. Most symmetric rank 2 distributions 

In this section we study the Tanaka algebra of a single Monge ODE 
and prove our main result. 

3.1. Calculation of the Carnot algebra. Let us consider equation 
(131) with fixed m < submanifold in the space of jets 

S C J™'"(R, M^) ~ M"+"^+=^(x, y,...,ym,z,zi,..., 

given as a hypersurface ym = F{x, y, . . . , ym-i, z, zi, . . . , Zn). Internally 
S ~ M"+'"+^(a;, y, . . . , Vm-i, z, zi, . . . , Zn) is canonically equipped with 
the rank 2 distribution 

C = {V^ = + vidy H h Fdy^_^ + zid^ H h Znd^^_^,d^J. 

Its (weak) derived flag has growth vector (2, 1, 2, ... ) if and only if 

In what follows we will assume this non- degeneracy condition which, 
according to Theorem [5l implies dimsym(C) < oo. 

In this case the derived distribution is C2 = C + {dz„_^ + Fz„dy^_^) 
and the second derived is 

= {T^x, dy^-i,d;,„, dz„_^,dz„_2 + Fz„dy^_2). 

The next (weak) derivecl^ are: 

^i+i — {^x, dym-11 • • • 1 dz„-, ■ ■ ■ dz„_i^^, dz„_i + Fz^dy^_.) , 

Cj+i = {V^, dy^_„ ...,dy, dz„, . . . dz„_^), i<m<j. 



These formulas imply that the weak and strong flags of C coincide, as was 
mentioned in Remark [TJ this plays a role in Section U) 
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Thus the dimension grows by two from grading —3 till grading — m — 2 
(resp. —m — 1 if = m) and then grows by one. The growth vector is: 




for m < n 



m n—m—1 



and 

(2, 1,2,..., 2,1) form = n. 

n-l 

Proposition 7. Carnot algebra of equation ^ does not depend on a 
point X E £ and is equal to the Lie algebra fm,n given by the following 
table of brackets (only upper- diagonal non-zero entries are shown). For 
m < n — 1: 





ei 




62 




64 64 . 


■ 6^+2 6m+3 (^m+A ■ ■ ■ 6„_|_i 


ei 




62 


63 


64 64 


65 6^ . 


■ 6m+3 6m+4 • • • 6^+1 












e', ■ 




62 















The GNLA structure is given by the rule: degCj = —i = dege^. 
For m = n the structure equations change to: 





6l 


e'l 


62 


63 


63 64 64 . 


6m 6m+l 6^_,_]^ 6^_(_2 


61 




62 


63 


64 


64 65 65 . 


■ 6m+l ^m+1 e^_^_2 


e'l 






e's 


e'. 


e'5 . 


■ ^m+l 6^^2 


62 















/Ind /or n = m + 1 t/ie table is: 





61 


e'l 


62 


63 e^j 


64 64 . 


■ 6m+l 6^^_j^ 6m+2 6^_^2 


61 




62 


63 


64 64 


65 e'g . 




e'l 










e'5 • 


e' 


62 















In other words, fm.n is obtained through m successive ID central 
extensions (this has sense to be called one central extension of order 
m) of the Heisenberg algebra 7in+2 = (ei, 6'^^, 62, 63, . . . , 6„+i). 

Proof. The generators are given by the formulae (we suppose F^^^n 7^ 0; 
we use the same symbol for the vector and its image under quotient to 
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GNLA) 

ei = -V^, e\ = e[ = F,^,^dy^_^^^ (2 < i < m + 2), 

ej = d,^_^^, + F,^dy^_^^, + c,9y,_^+2 (2 < J < m + 2), 

ek = d,„_,^, {m + 2<k<n + l) 

with coefficients C2 = 0, Cj = F,^_, + F,^Fy^^^_, - (j - 2)P,(F,J (3 < 
J < m + 2) playing no role in what follows. For n = m + 1 the last 
line in the above formula is void, and for n = m we shall additionally 
modify the range of the middle line (2 < j < m + 1). The structure 
equations follow. □ 

3.2. Calculation of the Tanaka algebra. Now let us describe the 
Tanaka algebra of equation ([3]) with Fz^z„ 7^ 0. 

Proposition 8. Prolongation tm,n = fm,n of the GNLA fm,n has the 
following graded structure (the number means the dimension while the 
subscript - the grading): 

(l(-m~2), 2(_m_i) . . . , 2(_4), 2(^.3), 1(_2), 2(_i), 2(0)) for n = m; 

then for n > m > 1: 

(l(-n-l) • • • l(-m-3)5 2(_m-2) • • • 2(-3), l(-2), 2(_i), 3(o), 1(1) . . . l(n-m-l)); 

and for n > 2, m = 1 we get: 

(l(-n-l), • • • , l(-4), 2(__3), 1(_2), 2(_i), 3(0), 2(1), 1(2), . . . , l(n-2))- 

Let us notice that for m = l,n = 2 the prolongation is different. 

Indeed, in this case fi,2 is isomorphic to the Lie algebra 02 with the 
grading 

(2(-3), l(-2),2(_i),4(o),2(i), 1(2), 2(3)). 
In the case m = n = 1 the prolongation is infinite, corresponding to 
the contact algebra of symmetries of J-'^(R, R): 

(1(_2), 2(_i),4(o), 6(1), 9(2), 12(3), 16(4), 20(5), 25(6), 30(7), 36(8) • • • )• 

To be more precise the GNLA structure of im,n = ®gi is given by 
the structure equations of Proposition [7] coupled with the following. 

Proposition 9. Forn = m: tm,m = fm,m®go> with the space go Abelian 
and generated by Cq^jCq^ given by 

[eo'',ei] = aei, [e'^%e[] = ce[, 

K,ek] = {{k - l)a + c) Ck, k>2, 

K,4] = ((^-2)a + 2c)e^, k > 2. 
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For n > m > 1 the space go C tm,n is a Borel suhalgehra of g\{g-i) 
and is generated by Cq^^, 6°^°, Cq^^ given by 

r abc 1 I L / r abc / 1 / 

[Cq ,ei\ = aci + be^, [Cq ,ej=cei, 

efc] = ((fc - l)a + c) efc + (fc - 2)6 e'„ A; > 2, 

K"',e',] = iik-2)a + 2c)el k > 2. 

Further on, tm,n = fm,n © S'o © " " " © fl'n-m-i; where gi for i > are one- 
dimensional and generated by C-i given by (only non-trivial relations 
are shown, and we omit commutators of non-negative gradation^): 

r 1 r 1 {fc-2)(fc-3){fc-4) / 
[e_2,eij = e_i, [e_2,efej = ^ '-^-^ e^_2, ... 

[Cm+I-Ji; Ci] = 6^+2-715 [^m+l-nj ^k] = (n-m) ^fc+m+l-n 

( in the above formulae for [e_p, e^] we restrict p + 3<k<m + p + 2). 

When n > 2, m = 1, the above formulae hold true with the only 
difference that gi is two-dimensional and is generated by ei.°,e°\, so 
that we change e_i to in the above, and add the only nontrivial 
relations 

[e°\, ei] = eo-^'°'^ = ef^ - 2 . [e°\, e\] = (1 - r^-) . 

[e°\,efc] = {n" -k'' + 2k-A) e^-i 

for the second generator of gi . 

Proof. Consider at first the case n = m. Then g_i canonically sphts as 
]Rei©Re'^, where ei is given (up to scahng) by the relation [ei, Cm+i] = 
(with Cm+i = ad™~"^(e2)) and e[ by the relation [€'^,63] = (with 
63 = ade'(e2)). Thus go is diagonal, and in fact it equals the space 
{bq'^} from the display formula of the Proposition. 

To see that gi = 0, consider its element e_i. It is uniquely given by 
4 parameters: 

[e_l,el] = eS^ [e_i,e;]=ef' . 
This implies the further relations: 

[e_i,e2] = ce'i - a'ei, [e_i, 63] = (a + 2c) 62, [e_i, 63] = (2a' + c') 62, 
[e_i, 64] = (3a + 3c) 63, [e_i, e^ = (2a' + c') 63 + {a + 2c) e'3, 

[6-1,65] = (6a + 4c) 64, . . . 
Relation = [6_i, [e[, e'^]] = (3a' + 3c') 6'3 yields c' = —a'. 



^which restore uniquely by the derivation property. 
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There are two ways to evaluate [6-1,65] since 65 = [61,64] = [6'i,64] 
and they do not match, implying a' = c' = 0, c = 0. We continue 
calculation and get 

r 1 (A;-l)(fe-2) r / 1 (A;-2)(fe-3) / 

[6-1, 6fe] = ^ — '-^ — -atk-i, [e_i,6fe] = ^ — '-^ — -ae^-i- 

For k — 171 + 2 (writing [6_i, em+i] = instead of em+2) we get a = 0. 

Let now n > m > 1. Then only the line Me'i C 51-1 is canonical, and 
it's easy to check that go is upper triangular in the basis 61,6^, with 
general element eg^'^ as in the display formula in the Proposition. 

Again to calculate gi we describe its element e_i via the formulae 

[6-1,61] = 6g^ [6-1,6;] = 60"''''^', 

which imply 

[6-1,62] = {c-b') e[ - a' 61, 

[e-i, 63] = (a + 2c - b') 62, [e_i, 63] = (2a' + c') 62, 
[6-1, 64] = (3a + 3c - b') 63 + 663, [e_i, 64] = (2a' + c') 63 + (a + 2c) 63, 

[e-i, 65] = (6a + Ac - b') 64 + 3664, ... 

Next = [6-1, [6^,63]] = (3a' + 3c') 63 yields c' = —a'. Evaluation of 
[6-1,65] in two ways using 65 = [61,64] = [6^,64] gives two different 
results, whence a' — d — 0, b' — c. 

The further calculations go without contradictions with the result 

[e-i, e.] = {i^^^a + {k - 2)c) 6,-1 + ^-^^^be',^,, 

[e-i,e',] = (M^a + 2(A:-3)c)e',_,. 

Now we need to specify and consider first n = m + 1. Substituting the 
last expression into [c-i, [61, 6^+2]] = and [6-1, [61, 6^_,_2]] = we get 
3 independent constraints implying a = 6 = c = 0. Thus in this case 
^1 = 0. 

Let now n > m + 1. Then the relations are [e_i, [ci, e„+i]] = and 
[6-1, [ei,6^_,_2]] = implying only a = c = 0, thus making C-i with 
6=1 the generator of gi. 

Next we consider g2 with element 6-2 given by 

[6-2,61] = a 6-1, [6-2,6'J = q;'6-i. 

This yields 

[6-2,62] = -Qi'e^^°, [6-2,63] = a'e[, [6-2,63] = 0, 

[6-2,64] = a' 62, [6-2, 64] =0, [6-2,65] =q;'63 + q;6'3, ... 
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Evaluation of [6-2,65] via 65 = [61,64] = [6'^, 64] gives la' and 0, 
whence a' = 0. Ifn>m + 2 there are no more contradictions and we 
get 

\p P 1 - (fc-2)(fc-3)(fc-4) / r ^/ 1 _ n 

[6_2,6fcJ — g tt6;j_2, [6_2,6^J— U, 

i.e. (72 is one-dimensional. 

For n = m + 2 the last relation in fm,n forces ct = and thus 5^2 = 
stops the prolongation. Continuing in the same way we uncover the 
result for n > m > 1. 

The case n > 2, m = 1 is considered in the same manner giving us 
the required prolongation fi Indeed, the 5^0 = {^o*'^} is the same with 
relations [8\ is the Kronecker delta) 

[6S^6l] = a6i + 66;, \e-^\(^^^ = ce^^^ 
\et\ ek\ = {{k - l)a + e) 6^ + 5l b e', {k>2), e',] = (a + 2e) e',. 
Next is Qi the elements of which are specified by the relations 

[6_i,6i]=6r, [6_i,6;]=6f'^'. 

Similar to the above we obtain a' = c' = 0, a = —2c, b' = {1 — n^)e. 
Thus the general element is 6^^ given by the only nontrivial relations: 

[e% ek] = {{k-l){3-k)+n'- l)c6fc„i + 6tbe',. 

Next comes g2 and the same arguments show that it is one-dimensional. 
Continuing we get Qp = (e-p) for 2 < p < n—2 with the only non-trivial 
relations (between elements of positive and negative gradings): 

[C-p, 61] = 6i_p, [6_p, 6p_|_3] = 63 

(and we mean e]^^ instead of 6_i in the very first formula for p = 2). 

Prolongation of fi^i is rather standard, and the case of fi^2 is due to 
Tanaka [T] (see also |Yj). □ 

Corollary 3. dim tm,n = 2n + 4 for n > m > 1 and it equals 2n + 5 
forn > 2,m = 1. In exceptional cases dimti^2 = 14 and dimti^i = 00. 

Thus in non-exceptional cases dim frn,n < 2 dim fm,n — 1 with the 
equality only for m = 1 (the deviation being 2m — 1 for m > 1). 

3.3. Symmetric models. In this section we shall calculate the inter- 
nal infinitesimal symmetry algebra sjm{Sm,n) for the model equation 
(jlj). The upper bound on the number of independent symmetries is 
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given by Corollary [3], and we shall show that this upper bound is re- 
alized by an explicit construction. An isomorphism between fm,n and 
sym{Sm,n) together with algebraic description of the latter will be given. 

If m > 1, the symmetry algebra has dimension dimsym(£^m,n) = 
2n + 4. These symmetries split into two types. The first type are the 
point symmetries with the generators: shifts 

^0 = d,, Yi = idyiO<t< m), Z, = (0 < J < n) 

and scalings 

Si = X dx + {m — 1) y dy + {n — z dz, R = y dy + \z dz- 

These form a subalgebra of symmetries of dimension m + n + 3, which 
acts transitively on £m,n (of dimension m + n + 2). As we shall see in 
Proposition [m the subalgebra {Sq, Si^Yi, Zj) C sym(£^m,n) coincides 
with the Carnot algebra fm,n, while the vector field Si + 2>R corresponds 
to the (Euler vectoi@ e G (7o- 

Then there are higher symmetries with generators (0 < < n — m): 

k 

^n+k — ^ ^1 \ p ) (k-p)\ ^"~™-P^J/ (n+fc)! 

p=0 

To verify that these are symmetries we prolong them to 

m k n 
7 — n \ ^ \ ^ ( -\\p (m+p~l-i\ x'^-P ^ a I \ ^ x"+''-^ a 

^n+k ^ ^) \ p ) ^n-m-p+i ^yi "T ^ {n+k-j)\ 

(with {^~^) = 5°) and check that 

Zn+kiVm - zl) = 0. 

When n = m, Zn is also a point symmetr}!^; otherwise Zn+k are genuine 
higher symmetries. 

Thus the totality of symmetries is 2n + 4 and since by Theorem H] 
dimsym(£^m < dimfm„ = 2n + 4, there are no more continuous 
symmetries. 

For m = 1, n > 2 there is one additional symmetr}|l, namely 
5*2 = x'^dx + {nzn-ifdy + {2n - l)xzdz., 
and then we hav^ dim sym(£^i^„) = 2{n + 3) — 1. Thus we have proved: 

^this is the unique vector characterized by the condition [e,v\ = fc • u V u G g^. 
'''this explains why this case separates in Proposition [T] 
^for n = 2 this symmetry persists as well, but there are 5 more symmetries. 
Vth dim£i,„ = n + 3, so that 2(n + 3) - 1 = dimS'r*Af, cf. [DZ] , 
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Proposition 10. The algebra sym(£^m^„) of internal symmetries of 
equation ^ has dimension 2n + A for n > m > 1 and dimension 
2n + 5 for n > 2,m = 1 with the basis Si, R, Yi, Zi. 

Let us describe the Lie algebra structure for sym(£^m,n)- Assume 
at first m > 1 and let tlie index / = ,2n — m. The non-zero 
commutators are the following (with Y_i = Z_i = 0). 

[5*0, 5*1] = 5*0, [5*0, Vj] = Fj-i, [5*0, 2';] = Zi^i, 

[Si, Y4 = {t + 1- m)Y„ [Si, Zi] = {l + \- n)Zi, 

[Y„R]=Yi, [Zi,R] = \Zi. 

These are accompanied by splitting Zi into Zj (0 < j < n) and Zn+k 
[0 < k < n — m) and the relations: 

[Zj, Zn+k] = 2{-l)'' {''"l~~^)Yj^n+m+k, U - ITL - k < j < U - k. 

For m = 1 (when Yi are reduced to solely Yq) we add the following 
commutators for ^2 

[^0, ^2] = 2Si, [Si, S2] = S2, [S2, Zi] = {{1 + 1- nY - n')Zi^i. 

By Theorem|l]and PropositionfTOlthe symmetry algebra of our model 
is isomorphic to the Tanaka algebra im,n = fm,n- We make this explicit. 

Proposition 11. The isomorphism sym(£^m,n) — tm,n is given by the 
following identifications: 

C _ c Q — 1 flOOl flioo Q — p — 1 oOOl 

•JO — Ci, 01 — 2% ~ ^0 ! ^2 — —G-i, -tt— — , 

Yi = {0<i <m), Zj = Cn+i-j {0 < j < n), Zn = -2e[, 

Zn+i = 2er, Zn+2 = -2e% Z^+k+i = 2(-l)Vfe {!< k < n-m). 

For m > 1 we shall remove S2 and e°\ with change ei^^ e_i. And 
form = n we identify e'^'^ ^ e'^ and remove Cg-^^ (Zn+i is non-existent 
in this case). 

In particular, we obtain the weights giving gradation of sym(£^m „): 
w{Sk) = k — 1, w{R) = 0, w{Yi) = i — m — 2, w{Zj) = j — n — 1. 

For m = 1, the radical is r = (Yq, Zi, R)o<i<2n-i and the vectors 
5*0, 5*1, 5*2 define an s/2 complement. The length of the derived series 
for r equals 3. The nilradical is {Yq, ^i)o<«<2n and its class (the length 
of lower central series) equals 2. 

For m > 1 the symmetry algebra is solvable. The length of the 
derived series equals 3 for n = m and 4 for n > m. The nilradical is 
{So, Yi, Zi)o<i<rnfi<i<2n-m and has class equal 2n-m + l. 
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3.4. End of the proof and a generalization. Now we can prove 
our main results. 

Proof of Theorem [T]. Let A be the distribution correspond- 
ing to equation (JSj) with non-degeneracy condition Fz„z„ 7^ satisfied 
and C the one for equation (jl]). We have demonstrated inequality 
dimsym(A) < dimsym(C). 

Moreover by virtue of Theorem H] the equality 

dimsym(A) = dimsym(C) = dimt^,! 

implies sym(A) ^ tm,n and the distribution A of the equation S is fiat. 

Therefore it is locally isomorphic to the standard model, which co- 
incides by the above with the distribution C. Thus the equation S is 
locally internally equivalent to the equation Sm,n ©• D 

Consider Monge ODEs with non-degeneracy condition F^^^n 7^ 
(sufiicient for dimsym(A) < oo). They correspond internally to rank 
2 distributions on a manifold of dim(£^) = m + n + 2 = k. Denote the 
class of all such distributions with m > fj, and fixed k by k). Our 
calculations imply: 

Corollary 4. Among distributions in the class D{m, k) the most sym- 
metric are locally equivalent to the distribution of equation Q). 

This implies in turn our Corollary [H Thus if we denote p^+s = fi,n, 
then the largest algebra of symmetries of distributions in S}(1, fc) is p^. 

Remarks: 1. Any rank 2 distribution can be realized as the Cartan 
distribution of an underdetermined ODE system. By the same trick as 
in D-modules this can be rewritten as a single higher order ODE. The 
equation manifold has however higher dimension and the former distri- 
bution is obtained from the latter by a sequence of de-prolongations. 
Thus without non- degeneracy condition Monge equations are internally 
indistinguishable from the general rank 2 distributions. 

2. Dimension 5 admits only one Lie algebra ps of the type (2, 1, 2). 
In dimension 6 there are already 3 GNLA of (2, 1, 2, l)-type: pg = 
fi,3! 1)6 = f2,25 and one more eltg not realizable as a distribution of a 
nondegenerate single Monge ODE |AKj . There is a reason to call them 
parabolic, hyperbolic and elliptic, whence our choice of the letters. 

3. In higher dimensions the number of GNLAs grows unboundedly, 
and it is impossible to trace the classification. Study of low-dimensional 
cases and some other types of underdetermined ODEs leads to the 
following claim, which will be discussed in a forthcoming paper |AK] : 
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Conjecture: Equation Si^k-s represents the most symmetric 
distribution C among all rank 2 distributions A on a manifold £ of 
dim(^^) = k with dimsym(A) < oo. 

4. Lie-Backlund type theorem 

Consider a bundle it : E ^ M. A Lie transformation of the jet space 
J'^TT is, by definition, a diffeomorphism preserving the canonical Cartan 
distribution C = {k is not the rank, but we write it as a superscript to 
avoid confusions). The classical Lie-Backlund theorem states that a Lie 
transformation E of the jet space J^ir is the lift of a Lie transformation 
/ of J'^TT, where e = 1 and / is a contact transformation of J^vr in the 
case rankvr = 1, while e = and / is a (local) diffeomorphism of J^vr 
(point transformation) for rankvi > 1. 

Various generalizations of this to Lie transformations of equations 
£ C J^TT are known |KLVt IAKO[ IGKj . In this section we present a 
generalization to the case of Monge equations ([3]) . 

4.1. Lie transformations of the mixed Jets. To prove our main 
result we first consider the case of mixed jeto 

jkiM(^M,Ni X A/'a) = J^'{M,Ni) J^'(M, A^a) 

with points {x , j^^ (si) , j^^ (32)) , where i'^^{si) denotes fcj-jet of the map 
Sj : M — i> Aj at the point x G M. There are natural projections 

TTi : J'^^''=^(M, A^i X N2) ^ J''(M, Ni) {i = 1, 2), 
^3 . jk^M^M.Ni X N2) j'''-^''"'-\M,Ni X N2). 

The Cartan distributions on the space of pure jets induce the one on 
the mixed jets by the formula 

Its (weak=strong) derived fiag is 

^(-fci.fe)^. _ (rfvr^)-i(C^i.'=2). 

Note that in the case considered in the next section, namely dim M = 
dim A^i = dim A^2 = 1, the growth vector ofC^^'''^ is (3, 2, . . . , 2, 1, . . . , 1). 

To exclude the known case J^{M,Ni x A'2) = J^'^{M,Ni x A'2) of 
pure jets, we assume in the next statement that ki < k2. 



^'^The results of this subsection have a straightforward generahzation to the case 
of more general mixed jets . 
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Theorem 12. 1. Every Lie transformation F of the mixed jet space 
j'=i''=2(M, Ni X N2) is the prolongation F^''^^ of a Lie transformation Fq 
of the space J^^^^-^^{M,Ni x N2). 

2. Under identification J^'^{M,Ni x N2) = x J^'(M,iV2) (for 
k = k2 — ki) every Lie transformation Fq is given by the formula 

where e = 1 and f : J^{M,N2) J^{M,N2) is a contact transforma- 
tion in the case dim Ai"2 = 1, while e = and f^''^ is the lift of a point 
transformation f : M x N2 ^ M x N2 for dim N2 > I. 

Proof. For 1 we mimic the proof of the classical Lie-Backlund theorem: 
from the derived flag of C''^ '^^ we deduce that the Cauchy characteristic 
space of {C'^'^'^^)j^i is Kerldir^) and thus every Lie transformation of 
jkiM f[iQQYS over a Lie transformation of J^^~i'^'^~j for any 1 < j < 
ki. Standard arguments show that the former transformation is the 
prolongation of the latter. 

To prove 2 let us notice that the Cauchy characteristics space of 
the Cartan distribution C°''^ equals TNi = Ker((i7r2). Thus every Lie 
transformation of J^'^[M, Ni x N2) fibers over a Lie transformation of 
the pure jet space J''{M,N2) and an application of the classical Lie- 
Backlund theorem finishes the proof. □ 

4.2. Lie transformations of Monge ODEs. Internal equivalence of 
two equations S,S' of type ([3]) is a diffeomorphism (p : S ^ S' which 
maps the induced Cartan distributions Cg to Cgi. Proposition [7| implies 
that two nondegenerate Monge ODEs with different pairs (m, n) are 
never internally equivalent. 

External equivalence of two Monge equations S,S' G J™''" (]R, M x R) 
is a Lie transformation of J™'"(R, RxM) mapping S to S'. In particular, 
it induces an internal equivalence. 

Theorem 13. Every internal equivalence of two non-degenerate Monge 
equations ^ with n > m > and {m,n) 7^ (1, 1), (1,2) is induced by 
a unique external Lie transformation. 

Contrary to other Lie-Backlund type theorems [KLVl [UK] . it is not 
enough here to study Cauchy characteristics of the derived distributions 
or characteristic subspaces. We shall use the Tanaka theory. 

Proof. Let C be the Cartan distribution of ([3]). Its (weak=strong) de- 
rived flag Ci was calculated in Section 13.11 

According to the proof of Proposition [9l the Cartan distribution C 
of (j3]) contains the canonical (from internal viewpoint) line C = {e'l). 



20 



IAN ANDERSON, BORIS KRUGLIKOV 



Namely the Carnot algebra fm,n contains precisely one vector f = e'^ G 
g^i (up to scaling) satisfying ad^ = 0. From external viewpoint this 
line C is generated by the vertical vector (in internal coordinates dz„, 
in external dz„ — F^^dy^). The quotient of £ by the corresponding 
foliation can be naturally identified with J'"~^'"~^(M, M^). Indeed, if 
we denote hj j : £ J"^'"^ the inclusion and hj n : £ ^ £/C the 
projection, then for : J"^'" jm~i,n-i have: n = o j. 

Furthermore, the distribution H = ("D^, dz„, <9z„_i, C C3 is also 

canonical (intrinsically defined) as it equals C2 © [e'i,C2] (in the Carnot 
algebra it corresponds to g_i © g_2 © (63)). In the quotient of £ by 
C the distribution H goes to the Cartan distribution of the jet space, 

d7T{n) = c'"-!'"-^ . 

Consider two Monge ODEs £,£' C J'"'"(M,M2) and an internal 
equivalence "if : £ ^ £'. It satisfies ^,(£) = £' and "^^(H) = W 
and thus induces a Lie transformation $ of the jet-space, i.e. we have 
the following commutative diagram: 



By Theorem [T2I $ equals the prolongation of a Lie transformation (j) of 
J°'"-'"(M, M^). To complete the proof we must show that ^! = 

First, we claim that maps £ to £'. Indeed, a point h E £ 

is represented by the class of tangent spaces T^(b)(7r(L^)), where 
runs through integral curves of {£,C) passing through h. Thus ^(6) is 
determined by the collection of spaces $*T^(b)(7r(L^)), which coincides 
with 

T^(.^,)0"^'\7:{L'))) = T^,(<^,„-.,(,))(7r'(</)("-)(Li))) 

due to equality $ = This implies that 0("~'")(£^) = £' , and 

let us denote 

J = ^-1 o :£ ^£. 

This map is a Lie transformation fibered over the identity map of 
jm-i,?i-ij^]^^ M^). Writing conservation of the Cartan form dzn-i — Zn dx 
we obtain that / is indeed the identity. □ 

Remark 2. The exact formulae from Section 1 5*. 31 show that symme- 
tries of 0) are induced by exterior transformations of jets J'"'" as the 
theorem claims. They have block form generalizing the feedback trans- 
formations in control theory (with variable z, control y and time x). 
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Let US discuss the exceptional cases. When m = n = 1, the in- 
ternal symmetry group of £ : yi = f{x,y,z,Zi) is isomorphic to the 
contact group of J^(M, M) and is parametrized by a single function in 
3 variables. On the other hand, the external symmetries can be cal- 
culated directly since, by the classical Lie-Backlund theorem, they are 
prolongations of point symmetries of J''(M, M^). In particular, for the 
model equation Si^i : y' = {z')^ the group of point symmetries is 10- 
dimensional, and thus is much smaller than the infinite-dimensional 
internal group. 

When m = 1, n = 2, the internal geometry of the equation S : yi = 
f{x, y, z, Zi, Z2) was described by E. Cartan |Clj . The internal symme- 
try group is at most 14-dimensional, and is realized by the equation 
£1,2 '■ y' = {z")^ (see |AKOt lAF] for the explicit generators). The 
symmetry group is the split real form of the exceptional Lie group 
G2. This group has a 9-dimensional subgroup of external symmetries 
coming from Lie transformations of J^'^(M, M^). But the remaining 5 
generators are not induced by the external symmetries. Thus, again in 
this case, the internal geometry is richer than the external geometry. 

5. INTEGRABLE EXTENSIONS 

In naive terms, a covering S or an integrable extension [KVl IBG] 
of a general partial differential equation S is given by introducing an 
additional dependent variable v and relations D{v) = F[x,u,v], where 
F are vector differential operators on the vector-functions u = u{x), 
V = v{x), such that the expanded system is compatible by virtue of 
S. In this case solutions of the new system S can be found from the 
solutions of S by ODE methods. 

In the case of (systems of) underdetermined ODEs with u = {y,z, . . .) 
the additional equation has the form 

v' = f{x,u,u\ ...,v) 

and no integrability constraints arise. 

In this section we prove that flat integrable extensions of fiat model 
equations are in bijective correspondence with central extensions of the 
corresponding Carnot algebras m. 

5.1. Central extensions of GNLAs. Any grading preserving auto- 
morphism of a fundamental GNLA is generated by its action on g^i, 
and therefore can be identified with a subgroup S) C GL(5f_i). Moti- 
vated by our interest in central extensions, we call it the gauge group 
of m. The Lie algebra i) = Lie{Sj) equals the subspace go <Z g = m of 
degree derivations of m. 
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A d- dimensional central extension m of a Lie algebra m is defined by 
the exact sequence 

^ a ^ m ^ m ^ 0, (5) 

where a is a d-dimensional Abelian algebra belonging to the center of 
m. We are interested in the case, when both m and m are fundamental 
GNLA, and a has a pure grading —k, i.e. a C g-k- 

For our applications the initial grading of GNLA is (2, 1, 2, ... ), and 
then we have 3 < A; < m + 1, where m is the length of m. An extension 
of maximal grading A; = m + 1 is 

m = 9-m~l © (g-rn © " " • © Q-l) = O © m 

(direct sum in the sense of vector spaces, but not Lie algebras) . 

It is well-known (see e.g. |SLt [F]) that the 1-dimensional non-trivial 
central extensions of a Lie algebra m are enumerated by the elements 
of the Lie algebra co homology if^(m). Zero element of the cohomology 
group corresponds to the trivial extension fn = m © R and we will 
exclude such cases, which for exact sequence dSJ means that it splits. 

The same proof shows that for a GNLA m central extensions of 
grading —k correspond to the graded part of the cohomology if|(m). 
More generally, rf-dimensional nontrivial extensions are given by d- 
dimensional subspaces in that is, elements of the Grassmanian 
GvdH^. Isomorphism classes of central extensions correspond to orbits 
of the gauge group action on Gr^ so that applying an automorphism 
to the cohomology class leads to an equivalent central extension. 

Thus we have proved 

Proposition 14. Let m be a GNLA with fundamental space g_i and 
the gauge group S^. Its d- dimensional nontrivial central extensions of 
grading k are in bijective correspondence with Gr^ i^|(m) /i^ . 

Let us notice that the space go is always non-zero, since contains 
the scaling corresponding to multiplication by on g^k- However this 
scaling acts trivially on the Grassman space Gr^ifKm). 

5.2. Integrable extensions of fiat models. For general distribu- 
tions A an integrable extension (or covering) is a submersion 

TT : (M, A) (M, A), such that 

(i) dx7^ : A ^ A is an isomorphism for any x G M, 

(ii) for the derived distribution we have A2 fl Kerrfvr = 0. 

Then by the Frobenius theorem, any integral manifold C (M, A) 
lifts to an integral manifold C (M, A) passing through a given point. 
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Indeed, the induced distribution A fl TS^v on Sat = vr N C M is 
Frobenius-integrable (and thus its integration is reduced to ODEs). 

Remark 3. To recast this into the terminology of |BG] let us choose 
a general connection on the bundle vr, i.e. a distribution of horizon- 
tal spaces H C TM such that A2 C H . This H defines a vertical 
Pfaffian system J C fl^{M) and the integrable extension is given by 
I = Ann(A), / = Ann(A), related by I = vr*/ + J. 

In order to relate integrable extensions to central extensions, we 
restrict to integrable extensions for which tangents to the fiber V of 
TT have pure grading —k in the Carnot algebra m. The integrability 
condition A2 fl Ker dir = is then equivalent to the requirement k > 2. 

If dim A = 2, then the integrability condition is vacuous, since we 
require the growth vector starts (2, 1, 2, ... ). In this case A; > 3. 

As usual we assume that the distribution A is completely non-holono- 
mic. An integrable extension vr will be called non-trivial if A enjoys the 
same property, which requires additionally only that iterated brackets 
span TV (and surely only in this case the above grading —k is defined). 

Theorem 15. 1. Let n : (M, A) (M, A) be a non-trivial integrable 
extension with d- dimensional fibers V . Let va denote the Carnot algebra 
at X E M and m the Carnot algebra at y = 7i{x) G M. Suppose that 
= TxV correspond to pure grading —k in m, k > 2. Then va is a 
d-dimensional central extension of m (in grading —k). 

2. Conversely, let m be a nontrivial central d-dimensional extension 
of a fundamental CNLA m in pure grading —k, k > 2, and let (M, A) 
and (M, A) be the corresponding flat model distributions. Then (M, A) 
is an integrable extension of{M,A). 

Remark 4. In other words there is a bijective correspondence between 
classes (modulo diffeomorphisms) of integrable flat extensions of flat 
distributions (M, A) and classes (modulo automorphisms of CNLA) of 
central extensions of m (nontrivial and of pure grading in both cases). 

Before we start proving this theorem, let us establish the functoriality 
of the Tanaka map A i-^ m. 

Proposition 16. Any smooth map vr : (M, A) (M, A) induces a 
morphism of CNLAs m^. tn^(a:) ■ 

Indeed, dn maps the derived fiags Aj Aj and hence 

dvr : g_i = Aj/Aj_i Ai/Aj_i = g_i. 

Since vr* = dir respects commutators, the map is a homomorphism. 
Now we can turn to the theorem about extensions. 
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Proof. 1. Due to Proposition [16] we have the following exact sequence 
of vector spaces 

^ t) ^ ffi ^ m ^ 0. 

Since vr* is an epimorphism of GNLAs it suffices to show that belongs 
to the center of m. 

Let u G Q-i and f G D. The assumption C g-k implies [u/u] G 
g-i-k- On the other hand, [u,v] equals [^,ri]x mod A^+i-i, where ^ G 
r(A;) is a vector field with value u at x and r] G T{V) is the vertical 
field with value v at x. Choosing a vr-projectible vector field ^ over some 
^ G r(Ai), we have [C,,r]] G r(V). This implies [u,v] G D n = 0. 

2. Consider a o?- dimensional central extension 

By the Tanaka theory the fiat distribution A can be locally identified 
with the standard distribution on the Lie group M corresponding to 
m, see |T]. Let (M, A) be the corresponding pair for fn. We take both 
M = M X V and M simply-connected, and as nilpotent groups they 
are topological vector spaces. 

The distributions A and A are spanned by the left-invariant vector 
fields from g^i and g^i respectively. Define the submersion vr : M — > M 
to be the group homomorphism, induced by the Lie algebras homomor- 
phism w. Then dn maps A onto A and A2 onto A2. Since k > 2, the 
distribution A2 is horizontal and hence tt is an integrable extension. 

2'. Let us give an alternative more constructive proof. 

Consider the integrable extension given by linearly independent co- 
homology classes [ai], . . . , [a^], where the closed 2-forms G A^m* 
have grading k > 2. These forms are not exact on m, but ti7*aj are 
exact on the extension m, whence zu*ai = dPi and is a basis of V*. 

This construction can be made explicit as follows. Realize a, as left- 
invariant 2-forms on M. Since M is contractible, at = dcTi G VL^{M). 
We can choose ai and linear coordinates Vi on V in such a way that 
the forms f3i = dvi — TT*ai are left invariant on M. 

Let A = Ann(/3i, . . . fl 7r^^(A). We claim that A2 does not 
intersect the fiber V. Indeed, choose u, w G r(A). Without restriction 
of generality we can assume that w cover sections v, w of A. Then 
ai{v,w) = thanks to the assumption k > 2, and so d[3i{w,v) = 
Pi{[v,w]) = 0. Since A2 is generated by the commutators [v,w], it is 
horizontal. Therefore tt is an integrable extension. 

Changing 2-forms within the cohomology class or applying an auto- 
morphism of m changes the distributions A to equivalent. □ 
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Theorem 17. Consider a flat underdetermined ODE S : F[x,u] = 
(u = u{x) is a vector function and F a vector differential operator) 
corresponding to a rank 2 distribution A with Carnot algebra m. Let 
xh be a nontrivial central extension of the GNLA m of grading k > 3. 
Then its integrable extension can be realized as an underdetermined 
ODE system (v has the same dimension as the extension) 

F[x,u] = 0, v' = G[x,u] 

In particular, if m = fm,n and S = £m,n is the model Monge equation, 
then the integrable extension can be realized as a pair of ODEs (which, 
in some cases, can be reduced to a higher order Monge equation) . 

Proof. To keep the notations simple we restrict to the case of a flat 
Monge equation. In Section [XT] we constructed coordinates {x,yi,Zj) 
on S (which are not the canonical coordinates on M = £^ considered as 
the nilpotent Lie group corresponding to m = fm,n) ■ On the integrable 
extension M let us use the coordinates adapted for the projection tt, 
that is the coordinates of M together with vertical coordinates vi. 

Let uji,uj[,uj2, ■ ■ ■ be the left-invariant 1-forms dual to the basis 
ei,e'i,e2, ■ ■ ■ on m chosen in Proposition [71 They satisfy the Maurer- 
Cartan equations duJi = 0, duj[ = 0, duj2 = uj[ A ui, ...Notice that 
ui = —dx,uj[ = dzn and A = Ann(a;j,a;^- : i,j >2). 

Thus we can express the forms defining the extension (as in the 
previous proof) via the co-basis: 

/3i = dvi - cnuji - civuj[ - C12UJ2 - . . . , 

where the coefficients Qo- are functions on M. By changing the 1- 
forms to cohomologous /3; 1— > — d[J ciy dZn] we can assume that the 
coefficient Qi' is zero. Thus Pi = dvi — gidx mod {uJi,Uj : i,j >2). 

This implies that the integral extension A = : i,j > 2) 

correspond to the underdetermined ODE system 

^M = (^(n))2^ vi = g,{x,y,z,z',...). 

For general ODE systems S the proof is absolutely similar. □ 

Returning to rank 2 distributions we observe that any fundamental 
GNLA m with gradings (2, 1, 2, . . . ) is a sequence of successive central 
extensions of the unique GNLA with gradings (2, 1, 2) and this implies 
Theorem [31 

5.3. Examples. In this subsection we demonstrate how Proposition 
[HI and Theorem [T71 work in practice. 
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5D^6D. Let us start with the unique GNLA fi^2 with the grading 
(2, 1, 2). Its muhiphcation table in the basis ei, e[, 62, 63, Cg is (we hst 
all non-zero commutators): 

[ei,e'i] = 62, [61,62] = 63, [61,62] = 63. 

This is equivalent to the structure equations in the dual co-basis: 

duji = 0, du'i = 0, du!2 = uj[ Aui, du^ = UJ2 /\uji, du'^ = UJ2 A uj[. (6) 

The group H'^{xn) = (of pure grading 4) is 3- dimensional, the gauge 
group is = GL2 = GL((7_i) and the quotient F{Hj)/S^ consists of 
3 points represented by the following 3 classes (called parabolic pe, 
hyperbolic f)6 and elliptic zUq 6D algebras, because if^(m) bears the 
natural conformal metric of Lorenzian signature (2,1), see |AK] ) . 

To calculate the integral extensions we need to realize the above 
abstract Maurer-Cartan forms. This is easy since they are dual to the 
vector fields 61 = —V^, e[ = 8^2, ^2 = d^^ + 2z2dy, 63 = d^, 63 = 2dy: 

Ui = —dx, uj[ = dz2, UJ2 = dzi — Z2dx, 

LJ3 = dz — Zidx, Lj'^ = \{dy — z^dx) — Z2{dzi — Z2dx). (7) 

pg: The co ho mo logy class oi a = uj^Auji corresponds to the extension 

duJi = uj-i AuJi (8) 
or, in the dual (vector) form, to the additional commutator relation 

[61,63] = 64. 

The flat model, calculated from (I7l) as A cji = —d{zdx), is the 1- 
dimensional integral extension obtained from Hilbert-Cartan equation 
by v' = z, whence the model Monge ODE is 

^1,3 : y' = wr- 

1)6 = f2,2- The cohomology class of a = A lj[ + u'^ A ui corresponds 
to the extension duj4^ = a or, in the dual form, to 

[61,63] = 64 = [61,63]. 

The fiat model, calculated from a.s a = djS with 

P = {zxZ2 - ^y)dx- Z2dx = -^ydx mod (u;2, ^^3, ^^3)' 

is the 1-dimensional integral extension obtained from Hilbert-Cartan 
equation hj v' = y (constant can be absorbed), whence the model 
Monge ODE 

£2,2 : v" = {z" f. 
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zUq. The cohomology class of a = u-s A ui + u'^ A u'l corresponds to 
the extension du^ = a or, in the dual form, to 

[61,63] = 64 = [61,63]. 

The flat model, calculated from (JTj) as a = with 

P = -{z+\zl)dx-\z2dy + \zldzi = -{z+\zl)dx mod (co'2, co'3, cuo), 

D 2 / D 

is the 1-dimensional integral extension obtained from the Hilbert-Cartan 
equation by v' = — (z + i^;!). The resulting systems is not representable 
as a single Monge and (after absorption of constants) equals 

y = [z ) , V = z - [z ) . 

6D-^7D. There are 8 different fundamental GNLA of dimension 7. 
We will consider only half of them - those which are extensions of pg. 
The structure equation of pe are (El) + ([8]). 

The cohomology group consists of two pure gradings H'^{Pq) = Hi © 
if| and both summands are 2-dimensional. The gauge group ^ is the 
(3-dimensional) Borel subgroup of GL(5f_i). Both quotients P(if|)/i3 
and P(if|)/i3 consist of two points and so there are 4 central extensions. 

(2,1,2,1,1)1. Consider at first central extension of maximal grading 
5. The fixed point of the action is given by the cohomology class 
[uj^ A uji\. The corresponding central extension is the GNLA py given 
by the additional commutation relation 

[61,64] = 65 

and the integral extension is 

y' = {z''?. 

(2, 1, 2, 1, 1)2. Another central 1-dimensional extension is obtained 
with the cohomology class [oj/^ A uj[ — uj^ A UJ2] (i^-orbit through it is 
open). The corresponding central extension is given by 

[6^,64] = 65 = [63,62] 
and the integral extension is given by the system 

y' = {z")\ v' = {z"f. 

(2, 1,2,2)1. Now consider the central extensions of grading 4. The 
fixed point of the ^-action is given here by the class [uj-^ Auj[ + uj'.^A uji\. 
The corresponding central extension is the GNLA f2,3 given by 

[ei,6^j] = 6^ = [6'i,63] 
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and the integral extension is given by the Monge equation 

y" = {z"'f. 

(2,1,2,2)2. Finally, another central extensions of grading 4 cor- 
responds to an open orbit of the i^-action and is given by the class 
[uj'2, A uj']\ . The corresponding central extension is given by 

[^n ^3] ~ ^4 

and the integral extension is given by the system 

y' = {z"r, v' = {z"'f. 

5.4. Extensions of parabolic models p„. The GNLAs pn = fi,n-3 
are of special interest as the most symmetric cases. It is obvious that 
is a central 1-dimensional extension of p„. We would like to study 
other possible central extensions of p„ of pure grading. 

Notice that any fundamental GNLA with dimg^.i = 2 can be ob- 
tained from (2, 1, 2) by successive d-dimcnsional central extensions of 
only maximal grading. In studying the latter for the sequence of par- 
abolic models we observe the following Z2-periodic pattern. 

Theorem 18. The Lie algebra p„+i is a maximal grading central ex- 
tension of pn- For odd n such an extension is unique, while for even 
n — 21 there are two maximal grading extensions pn+i and p'^+i- 

The second is however less symmetric (meaning dim p'^_^_i < dimp„+ij 
and has no maximal grading central extensions. 

Proof. The Maurer-Cartan equations of p„ are (the index denotes the 
grading): 

duJi = 0, du[ = 0, dLJ2 = u[ A Ui, doj^ = UJ2 /\ oJi, du'^ = 002 Au[, 
dui = a;3 A cui , duj^ = a;4 A tui , ... du!n-2 = '^n-s A cui . 

The 2-cocycles in the maximal grading are 

-^n-l = A Ui, UJn-2 A Uj[ - UJn-Z AUJ2+ OUn-4 ACU3 ± O;/ A CUi-i) 

for even n ~ 21 and Zl_^ = {uJn-2 A iOi) for odd n — 21 + 1 

There are no coboundaries in this grading and so T'H^_^ — IRP^, 
where e = for odd n and e = 1 for even n. The action of 3-dimensional 
gauge group has two orbits in the latter case: fixed point [wn-2 Acui] 
and the complement, so that we can take the the second form of Z^_^ 
for n = 2Z as a representative. 



RANK 2 DISTRIBUTIONS OF MONGE EQUATIONS 



29 



Extension through the first class is possible for any n and the result is 
p„+i. Extension via the second class yields the GNLA p^+i for n = 21 
with the following additional commutation relation: 

[e[, en-2] = Cn-i, [e2, e„_3] = -e„_i, [es, e„_4] = e„_i . . . 

It's easy to check that there are no 2-cocycles Z2i{p'2i+i) = O5 so 
-^2«(p2/+i) ~ 0- Thus no further maximal central extensions exist. 

Finally to calculate the Tanaka prolongation of P2/4.1 notice that both 
directions Rei and Re'^ are canonical (invariant under automorphisms) 
in g_i because they can be uniquely characterized by the properties 
adei \g2_21 = 0, ade'^ \g_^ = 0. Therefore the space go is 2-dimensional 
Abelian (diagonal 2x2 matrices). Calculations similar to that of the 
proof of Proposition [9] show that g^i = and so dimp2/+i = 2/ + 3 is 
smaller than dim P21+1 = 4/ — 1 . □ 

The fiat ODE for p„ is given by the model equation Si^n-s- The 
flat ODE for is not realized as a Monge equation, but only as a 
system. To find it one follows the scheme of Section 15.31 and therefore 
obtains the following (non-standard) integral extension of 81^21-3'- 

6. Applications to integration of PDEs 

We conclude with some examples which illustrate how our study of 
rank 2 distributions for Monge equations can be applied to the inte- 
gration of overdetermined and determined PDEs. 

6.1. Prom underdetermined ODEs to overdetermined PDEs. 

The examples here generalize those of Cartan |Clj . Specifically, Cartan 
considered compatible systems consisting of pairs of scalar 2nd order 
PDEs on the plane with one common characteristic (which we denote 
as 2ii^2)- The equation S C J^(M^,IR) from the internal viewpoint is a 
6-dimensional manifold Nq equipped with a rank 3 distribution. This 
distribution always has a Cauchy characteristic and can therefore be 
reduced to a rank 2 distribution on a 5-dimensional manifold N^,, which 
is necessarily equivalent to a single Monge equation 

Now we consider the case of mixed order equations S consisting of 
pairs of compatible 2nd and 3rd order PDEs {E2 + E3). Again it can 
be seen that such systems can be reduced to rank 2 distributions, now 
on 6-dimensional manifolds, but we would like to emphasize that the 
reduction is slightly more complicated than in Cartan's case. 
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Indeed the system 8 is represented internally as an 8-dimensional 
manifold Mg equipped with a rank 3 distribution Q. This distribu- 
tion has one Cauchy characteristic, while its derived distribution Q2 
has two. Thus the pair (Mg, Q) has two reductions - one by Cauchy 
characteristic for Q (quotient) and one by Cauchy characteristic for 
the derived distribution Q2 (de-prolongation). 

These two reductions commute, so that we can do them in any order, 
and in Examples 1 and 2 we perform both at once. But for the pur- 
poses of Example 3 we shall need first to de-prolong Mg to M7 (with 
distribution Q of rank 3), and then to quotient (My, Q) to (Mg, A) 
(rank of A is 2). These steps are summarized in the following sequence 
of maps: 

(Mg,Q)-^(M7,Q)^(M6,A). 

In the first two examples (Mg, A) turns out to be one of the 6- 
dimensional symmetric models from Theorem [1] (this can be detected 
a priori via calculation of the symmetry algebra and concluding that 
the distribution is fiat in the Tanaka sense). Then we use their integral 
curves to integrate the corresponding PDE systems. 

Example 1: E2 + E3 ^ £^1 3. Consider the PDE system, written 
parametrically as 

"^xx T X , T X , '^xyy "^yyy ~^ 5'^ ^' (^) 

It is compatible and has a common characteristic ^ = — ^rVy. 

Internally the equation is Mg = M.^{x,y,u,Ux,Uy,T,Uyy,Uyyy) and 
the derived distribution Q2 has a two-dimensional space of Cauchy 
characteristics 

n={i = dx- \Tdy + {Ux - \TUy) du - \t''x^ dux 

+ K4rV - ^Uyy) du, + \T^Xdu,, - ^9., du^J. 

The distribution 11 is integrable and the quotient by its leaves 

Me = Mg/n ~ n\t, V, w, wu W2, Ws) 
can be identified with the space of invariants 

f — It -\- —n~ 'y ?i — 256 ^ , i 128 ^5 ^3 
b y ^ ^1 d., u 625 ^ ~'~ 625 ' 

w = — Ir^x"^ + yT^a;^ Uyy + x Uy — x Ux + u, (10) 

Wl = Uy + ^TXUyy — 3T X , W2 = Uyy ^ X , W3 = — X. 

The quotient distribution on the manifold Mg equals 

A = {dt + wldu + widu, + W2du,^ + w^dui^.d^,^) 
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and this corresponds to the model Monge equation £^1^3: 

v' = {w"'f. 

The inverse image of the general solution to this equation with re- 
spect to the projection map vr : Mg — » Mg given by (fTOj) (this introduces 
only 1 more parameter: vr"^ decomposes into preimage and prolonga- 
tion) yields the general solution of ([9]), written parametrically as 

^ = -(|)V'w, y = t-sr{t), 
u = m - ^ntr + jnt) ntf - s nt) nt) - '^nt) g{t\ 

where g'{t) = f"{tf. 

Example 2: E2 + E3 £2,2- Now let's consider the system 

'^'^xx y ) '^^yy ' ^xyy ^yyy ^ ' ('^■^) 

which is again compatible with common characteristic ^ = V^. — r'^'Dy. 

Internally, the equation is Mg = M.^{x,y,u,Ux,Uy,T,Uxy,Uyyy) and 
the derived distribution Q2 has a two-dimensional space of Cauchy 
characteristics 

Il={i = d,- dy + K - Uyr^)du + - u.,yr^)du^ 

+ (u^y - k)duy + iy- r)du,^ - T^dr, du^^y). 

The distribution 11 is integrable and the quotient by its leaves 

Me = Mg/n ~ R\t, V, vi, w, wi, W2) 

is identified with the space of invariants 

t = y — T, f = lnr — l + 2n — 2TUy + 2uxy — 2xy — r x^y + 2r x+ 

-I — — — 2x Ux + 2t X Uxy, vi = xy + - — Tx — 2x Uxy + 2uy, 

w = ^ - y - T xy + r + ^ + T Uxy - Ux, 
Wi = xy — T X — Uxy, W2 = —7 + X. 

In these coordinates the quotient distribution A on Mg coincides with 
the rank 2 distribution of the model Monge equation 82,2'- 

II I II \2 

V = (w ) . 
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Again the general solution of this equation provides the general solution 
to system (ITTl) as 

2s s 2s 2 

where g"{t) = fitf. 

We remark that the 2nd order PDE in ([9]) x^u^^ = u^y is a new 
example of a Darboux integrable equation. On the other hand, the 
2nd order PDE in (fTTj) 2Ma;a; + {2uyy)~^ = y"^ is not Darboux integrable 
and can therefore only be solved by compatibility technique. 

Example 3: E2 + E3 ^ 2E2. In this example we show that our 
theory of integrable extensions for Monge equations can be applied to 
construct integrable extensions of overdetermined PDE systems. In- 
formally speaking, if £ and £' are two PDE systems whose reduction 
by Cauchy characteristics and de-prolongations are rank 2 distribu- 
tions corresponding to equations £ and £' , then an integrable extension 
if : £ ^ £' can be lifted to an integrable extension : £ ^ £' . 

To be more precise, in the notations adopted at the beginning of 
this section we have the following commutative diagram (we mark by 
primes the canonical distributions of the manifolds Ni) 

(M7, Q) {Ne, Q) 



(Me, A) (iV5,A') 

where the horizontal arrows are the integrable extensions and the ver- 
tical arrows are the quotients of the rank 3 distributions Q and Q! by 
their Cauchy characteristics. 

We illustrate this construction with the example of Cartan's involu- 
tive symmetric model [Clj 

Uxx = |r^ UxY = \r\ Uyy = r. (12) 

This system reduces to ^1,2: V = {W"f with V = V{T), W = W{T) 
and the projection vr' along the Cauchy characteristic defined by 

T = -Uyy, V = -2U+2XUx+2YUy-Y'^Uyy-XYU^y-IX^U^y^ 
W = Ux- UyUyy + \YU^y + 
Wi = Uy- YUyy - \XU^Y^ W2 = Y + XUyy 
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Basing on the fact that £^1,3 is an integrable extension of 81^2 and on 
the reduction from Example 1 of this section, we will show that ^ is 
an integrable extension of (fT2l) . 

Recall that (M7, Q) is obtained from (Mg, Q) by de-prolongation via 
a Cauchy characteristic for Q2 not belonging to Q. Thus Q still carries 
a Cauchy characteristic and we get (M7, Q) ~ (^1^3, C^'^) x (R, R) and 
similarly (A^g, Q') — {£1 2, C^'^) x (M, M) for the Cartan example, because 
(Me, A) ~ (^i,3,Ci'3) and {N,,A') ^ {S^,2,C''')■ 

The central extension w : pQ ^ induces the integrable extension 
= X 1 : (£1,3, C) X (M,M) ^ (^1,2, C) X (M,M), where 

V9(t, f , W, Wi, W2, W3) = (t, V, Wi,W2, Ws) = (T, V, W, Wi, W2). 

By this we construct the required map $ in the diagram, and then 
we compute its inverse <l>~^ given by the formulae 

X = -{IX)\y + X Uyy)-\ y = -Uyy - |X, 
= {l)\-lX'U^y - \X{2Y + X^) UIy - {X^ + Uyy 

+2X Ux + 2YUy-2U- |Xr2) , 
Uy = |X U^y + (f + lY)U^Y + (i^' -Uy + \XY)Uyy 

+Ux-IXUy + ^XX 

where (X, Y, U, Ux, Uy, ■ ■ ■) are jet-coordinates on the equation Nq (IT^ 
and (x, y, u, u^, Uy, . . . ) are jet-coordinates on the equation M7. Pro- 
longation of the latter yields the required integrable extension Mg. 

Notice that U is not expressed in the above formulae but is obtained 
via quadrature, as is always the case with flat integrable extensions. We 
leave it to the reader to check that the integrability condition Dx{uy) = 
Dyiux) holds by virtue of f|T2l) and satisfies ([9]). 

6.2. Relation to the method of Darboux. According to [AFV] any 
Darboux integrable single scalar PDE on the plane £ C J^(R^,R) is 
a group quotient of a pair of rank 2 distributions with common freely 
acting symmetry group G, that is 

(^,C) = (Mi,Ai)x(M2,A2)/G. 

Notice that in this case dim£^ = 7, dimC = 4, but dimMj can be 
sufficiently large, depending on the order of intermediate integrals of £. 

Our theory of rank 2 distributions of maximal symmetry restricts 
the possible types of the factors. 

Theorem 19. Suppose that {Mi, Ai) have Carnot algebras fmi.nj- Then 
the pairs (mi, ni), (m2, ^^2) with normalization < rrij < Ui, n2 — ni = 
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k > 0, and nii < m2 for k = satisfy the condition 

mi + 1712 + k < 7 + 6, 

where 6 = 1 if mi = m2 = 1 or mi = 1, k > 0, and 6 = else. Thus 
there are no more than 53 different admissible triples {mi,m2, k) . 

Proof. The equality 

dim Ml + dim M2 — dim G = dim S = 7 

implies dim G = rii + mi + n2 + m2 — 3. Since G acts as the symmetry 
group for both (Mj, Aj) and since the symmetry group of a distribution 

is bounded from above by Tanaka algebra t^.n = fm,n, i-e. dimG < 
dimtm-^„- we get that minjdim tm-_„- : i = 1, 2} > ni + mi+n2 + m2 — 3. 
From Section 13.41 we know that dimtm^ n^ = 2nj + 5 for mj = 1 and 
2ni + 4 for > 1. The claim follows. □ 

Recall that distributions Aj corresponding to nondegenerate Monge 
equations ([3]) have Carnot algebras fm^.n^ The estimate of the theorem 
is rough and can be sharpened by a detailed analysis of the admissible 
(Mj, Aj) which actually arise for various types of equations £. 

Example ([AF]). The Goursat equation u^y = a{x, y)^UxUy is Dar- 
boux integrable if and only if both factors (Mj, Aj) can be realized as 
rank 2 distributions of the Monge equation 

y' = F(x,^,^',...,z(")) 

and the group G is either Abelian or 1-step solvable of dimension 2n+l. 
In particular, for the Darboux integrable equations 

2n 

both factors coincide with (the prolongation of) the model equations 
£i^n- In this case, the Lie algebra of G is equal to the nilradical r of 
the Lie algebra ti,„ given in Section 13.31 
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